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Abstract 

Spina bifida is one of the most common neural tube defects (NTDs) with a complex etiology. Variants in planar cell polarity 
(PCP) genes have been associated with NTDs including spina bifida in both animal models and human cohorts. In this study, 
we sequenced all exons of CELSR1 in 192 spina bifida patients from a California population to determine the contribution of 
CELSR1 mutations in the studied population. Novel and rare variants identified in these patients were subsequently 
genotyped in 190 ethnically matched control individuals. Six missense mutations not found in controls were predicted to be 
deleterious by both SIFT and PolyPhen. Two TG dinucleotide repeat variants were individually detected in 2 spina bifida 
patients but not detected in controls. In vitro functional analysis showed that the two TG dinucleotide repeat variants not 
only changed subcellular localization of the CELSR1 protein, but also impaired the physical association between CELSR1 and 
VANGL2, and thus diminished the ability to recruit VANGL2 for cell-cell contact. In total, 3% of our spina bifida patients carry 
deleterious or predicted to be deleterious CELSR1 mutations. Our findings suggest that CELSR1 mutations contribute to the 
risk of spina bifida in a cohort of spina bifida patients from California. 
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Introduction 

Neural tube defects (NTDs) are among the most severe and 
common of all human birth defects. The most frequent types of 
NTDs are spina bifida and anencephaly. The etiology of NTDs is 
complex and involves both environmental and genetic factors. 
Periconceptional folic acid supplementation reduces 50% to 70% 
of newborn NTDs [1]; however, the mechanisms underlying this 
protective effect remain unclear. In terms of genetic underpin- 
nings, monozygotic twinning and single gene disorders have long 
been associated with increased risks of NTDs [2]. Numerous 
exploratory candidate gene studies have highlighted a variety of 
biological pathways such as the folate and one carbon metabolism 
and transport [3] , DNA repair [4] , retinoic acid receptors [5] , and 
the planar cell polarity (PCP) signaling network [6]. 

The PCP pathway controls the polarity of cells within the plane 
of epithelium in both vertebrates and invertebrates. The PCP 
genes, including frizzled, dishevelled, vangl, flamingo (Celsr) , prickle and 
diego, were initially identified in Drosophila, and are highly 
conserved throughout evolution [7]. PCP signaling is required 
for the initiation of neural tube closure in higher vertebrates [8] . In 
mice, mutations in Vangl2, Celsrl, Dishevelled and Frizzled result in 
the NTD known as craniorachischisis [9] . In humans, mutations in 
FRIZZLED6 [10], DISHEVELLED [11], VANGL [12,13] SCRIB 
[14] and PRICKLE1 [15] have been associated with NTDs. 

CELSR 1 is well known for regulating the establishment and 
maintenance of planar cell polarity. During mitosis, CELSR 1 
recruits VANGL2 and FZD6 to endosomes. Following mitosis, 
CELSR 1, VANGL2 and FZD6 are recycled to the cell surface to 



re-establish cell polarity [16]. In Drosophila, CELSR 1 mediates 
homotypic interaction between adjacent cells and transmits 
instructive PCP signals. In zebrafish, the knocking down of Celsrl 
produced convergent extension (CE) defects [17]. In mice, Celsrl 
mutants exhibited craniorachischisis, the most severe form of 
NTDs [18]. In humans, functional CELSR1 single nucletotide 
variants (SNVs) have been identified in fetuses with craniorachis- 
chisis, and predicted-to-be-deleterious CELSR1 SNVs have been 
detected in a few cases with NTDs or caudal agenesis [19,20]. 
However, the contribution of CELSR1 mutations in the etiology of 
spina bifida is still unknown. In this study, we investigated the 
CELSR1 coding region sequence among a cohort of spina bifida 
infants born in California by Sanger sequencing. Further, we 
conducted in vitro functional analyses to validate the functional 
effect of the identified mutations. 

Materials and Methods 

Ethics statement 

The approval process includes detailed review by the State of 
California Committee for the Protection of Human Subjects (the 
primary IRB). All samples were obtained with approval from the 
State of California Health and Welfare Agency Committee for the 
Protection of Human Subjects. 

Human subjects 

Data were obtained from a population-based case-control study 
conducted by the California Birth Defects Monitoring Program 
(CBDMP). The CBDMP is an active, population-based surveil- 
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Figure 1. TG dinucleotide repeats variants in spina bifida. A: Sequence trace of control (top) and C.5719-5720delTG (bottom). B: Sequence 
trace of control (top) and TG duplication (bottom). C: Schematic representation of the CELSR1 predicted protein structure (accession number 
Q9NYQ6) with the domains and approximate position of TG repeats variants. 
doi:1 0.1 371 /journal.pone.0092207.g001 



lance system for collecting information on infants and fetuses with 
congenital malformations, which has been described elsewhere 
[21]. Included for study were 1 92 infants with isolated spina bifida 
and without other major birth defects (cases) and 190 non- 
malformed infants (controls). Cases randomly selected from all live 
born infants with spina bifida and a random sample of non- 
malformed infants were ascertained by the CBDMP correspond- 
ing to birth years 1983-1999. Among the 192 spina bifida cases, 
82 are White NonHispanic, 54 are native US born Hispanics and 
56 are foreign born Hispanics. Among the 190 controls, 81 are 
White NonHispanic, 54 are US born Hispanics and 55 are foreign 
born Hispanics. All of the 192 spina bifida included in this study 
are cases of myelomeningocele. The case and control infants were 
linked to their newborn screening bloodspots, which served as the 
source of the gDNA used in these studies. Bloodspots are collected 
on all newborns in California for genetic testing purposes by the 
State of California. The State retains the residual, unused portion 
of the bloodspot and makes these bloodspots available to approved 
researchers. 



DNA sequencing 

Genomic DNA was extracted using the Puregene DNA 
Extraction Kit (Qiagen, Valencia, CA) and amplified using the 
GenomiPhi Kit (GE Healthcare). Coding exons and flanking 
exon-intron regions of the human CELSR1 gene (NM_0 14246) 
were amplified by polymerase chain reactions (PCR) from the 
whole genome amplification (WGA) product. Primer sequences 
are available upon request. The PCR products were sequenced 
using the Prism Bigdye Terminator Kit (v3) on an ABI 3730XL 
DNA analyzer (Life Technologies, Carlsbad, CA). Both case and 
control samples were sequenced with either a specific forward or 
reverse primer. Sequencing results were analyzed using the 
Mutation Surveyor software V4.0.7 (Softgenetics, Stage College, 
PA). Detected mutations were subsequently confirmed by a second 
round of whole genome amplification, PCR and sequencing 
analysis. GenomiPhi Kit takes advantage of Phi29 DNA 
polymerase, which produces high fidelity during DNA replication 
due to its proofreading 3 '-5' exonuclease activity. The reported 
error rate of Phi29 is between 3xlO" B [22] to 5xl0" 6 [23]. The 
probability to generate the same artifact mutation in two rounds of 
WGA is 9 x 1 0 " 1 2 to 2 . 5 x 1 0 " 1 1 . In our CELSR 1 mutation screen, 
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A 




GFP-CelsrlWT GFP-Celsrl TGins GFP-Celsrl TGdel 

B 

GFP-Celsrl-WT GFP-Celsrl-TGdel GFP-Celsrl-TGins 



Figure 2. Subcellular localization of GFP-Celsrl wild type and TG repeat variants. A: MDCK II cells were transfected with GFP-Celsrl 
plasmids. Each image shown is representative of at least 50 examples. It demonstrated that the C.571 9-5720delTG and the C.5050-5051 insTG disrupt 
GFP- Celsrl membrane localization. Scale bar, 25 |rm. B: Western blot of GFP-Celsrl wild type and the two indels mutants. It showed that the c.571 9- 
5720delTG and the c.5050-5051 insTG changed the size of GFP-Celsr1 protein. 
doi:1 0.1 371 /journal.pone.0092207.g002 



we sequenced 192 spina bifida cases with 9045 nucleotides in the, 
CELSR1(NM_0 14246) gene. In total, we screened 1.74 xlO 6 base 
pairs. The probability to detect the same coding region artifact 
mutation in our study is less than 4.35 x 10 . 

Plasmids 

Mouse Ceslrl cDNA cloned into a pEGFP-Nl plasmid 
(pEGFPNl -Celsrl) was kindly provided to us by Dr. Elaine Fuchs 
(The Rockefeller University, New York, USA). Celsrl open reading 
frame (ORF) was sub-cloned to pEGFP-Cl. Human influenza 
hemagglutinin (HA) tagged VANLG2 (HA-VANGL2) plasmid was 
obtained from Dr. Hongyan Wang (Fudan University, Shanghai, 
China). VANGL2 ORF was sub-cloned into the pDs-RedCl at 
Xhol and Sail restriction sites. CELSR1 nonsense and missense 
changes were introduced into pEGFPNl -Celsrl by QuikChange 
II Site-Directed Mutagenesis Kits (Agilent Technologies, In- 
c.CA,USA). All plasmids were validated by sequencing analyses. 

Subcellular localization 

MDCK II cells were purchased from Sigma-Aldrich and 
cultured according to the manufacturer's protocols. One day 
before transfection, cells were seeded in 4 chamber 35 mm glass 
bottom dishes (4x 104/chamber) (In Vitro Scientific, Sunnyvale, 
CA). Plasmids transfection was performed using GeneTran III 
Tranfection Reagent (Biomiga, San Diego, CA) according to the 
manufacturer's manual. Forty eight hours later, cells were washed 
twice with PBS and incubated 5 minutes with Hochest 3342 



(1 ug/ml) (Invitrogen), then washed 3 times with PBS and fixed in 
4% PFA (paraformaldehyde in phosphate-buffered saline) for 
10 minutes at 37°C, followed by 3 times PBS wash. Cells were 
examined and photographed by an LSM710 laser scanning 
confocal microscope (Leica). 

Immunoprecipitation and immunoblotting 

HEK293T cells were grown and maintained in DMEM 
supplemented with 10% fetal bovine serum (FBS) on a 6-well 
plate at a concentration of 6 x 1 0' cell/ well before the day of 
transfection. 2 |Ag of GFP-Celsrl or its related mutant plasmids 
were co-transfected with 2 |Ig of HA-Vangl2 by Lipofectamine 
reagent (Invitrogen). Twenty-four hours post-transfection, cells 
were washed twice with ice-cold PBS and lysed with 300 ul lysis 
buffer. Lysate was pretreated with protein A/ G agarose, and then 
immunoprecipitated with 1-2 H-g anti-HA antibody and protein 
A/ G agarose at 4°C overnight. After washing three times with lysis 
buffer, the precipitates were run on SDS-PAGE followed by 
Western blot detection immunoblotting with the anti-GFP 
antibody. 

Results 

Novel rare mutations identified in CELSR1 among spina 
bifida patients 

The human CELSR1 gene coding region sequence contains 46 
TGTG and 3 TGTGTG dinucleotide repeats (Figure SI). Two 
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Figure 3. Subcellular localization of Celsrl and VANGL2 in MDCK II cells. A: The cells were transfected with GFP-Celsrl or DsRed-VANGL2. 
GFP-Celsr1 itself was localized to cell membrane while DsRed-VANGL2 was unable to localize to cell membrane. B: The cells were co-transfected with 
dsRed-VANGL2 (Red) and GFP-Celsr1 wild type, C.5719-5720delTG, C.5050-5051 insTG. In the presence of Celsrl WT-GFP, DsRed-Vangl2 changed its 
localization to cell-cell borders only when two Celsrl -expressing cells are in direct contact. Each image shown is representative of at least 30 
examples. It demonstrated that CELSR1 C.5050-5051 insTG and C.5719-5720delTG mutants failed to recruit DsRed-VANGL2 to cell-cell contact. 
doi:1 0.1 371 /journal.pone.0092207.g003 



TG dinucleotide repeat variants were identified in spina bifida 
cases (N= 192), and both of them were absent among the 190 
control samples (Table 1 and Figure S2). One was a TG-insertion 
(c. 5050-505 linsTG) and one was a TG-deletion (c.5719- 
5720delTG), both of which created a stop codon in the middle 
of the CELSR1 ORF. The C.5050-5051insTG created a stop 
codon at the 1706 th amino acid, and the C.5719-5720delTG 
produced a stop codon at the 1944 th amino acid (Figure 1). We 
also identified 1 1 missense SNVs in NTD samples but not in any 
controls, six of which were predicted to be deleterious or damaging 
by both SIFT and PolyPhen (Table 1). The six mutations are 
p.Arg2497Cys, p.Arg2354Cys, p.Glyl410Arg, p.Thrl362Met, 
p.Ilell24Met and p.Alal023Gly. Four of the 6 SNVs were 
predicted to be damaging and disease causing by MutationTaster 
and FATHMM, they were p.Arg2497Cys, p.Thrl362Met, 
p.Ilel 124Met and p.Alal023Gly. These four mutations were 
mapped to different domains of CELSR1: p.Alal023Gly and 
p.Ilel 124 were mapped to the cadherin repeats, p.Thrl362Met 
was mapped to a EGF-like domain and p.Arg2497Cye was 
mapped to the transmembrane domain of CELSR1 (Table 1). 
Among these four mutations, one (p.Arg2497Cys) was identified 
once in the 1000 genome sequencing project. None of the 11 
SNVs were previously detected in the Exome Variants Project 
(Table 1) or the two previously published NTDs' CELSR1 
mutation screen studies [19,20]. All rare mutations identified in 



this study have been uploaded to LOVD website (http:/ /www. 
Iovd.n1/3.0/home). 

The CELSR 7 C.5050-5051 insTG and the C.5719- 
5720delTG disrupt CELSR1 membrane localization 

Several point mutations (p.Ala773Val, p.Arg2438Gln, 
p,Ser2964Leu and p.Pro2983Ala) identified in humans craniora- 
chischisis cases were found to alter membrane localization [19]. In 
our study, both the C. 5050-505 linsTG and the C.5719- 
5720delTG variantsintroduced stop codons ahead of the CELSR 1 
transmembrane domain. Therefore, we predicted that these two 
variants would also impair CELSR1 membrane localization. In 
the CELSR1 subcellular localization assay, both the C.5050- 
5051insTG and the C.5719-5720delTG misplaced CELSR1 in 
MDCK II cells. Unlike wild type CELSR 1 which localizes almost 
exclusively on the cell membrane, CELSR1 protein with C.5050- 
5051insTG or C.5719-5720delTG distributed throughout the 
cells, including the cytoplasm and the nucleus (Figure 2) in all of 
the cells examined (100%). 

The CELSR1 C.5050-5051 insTG and the C.5719- 
5720delTG disrupt VANGL2 cell-contacts localization 

It was demonstrated that Celsrl is required to recruit Vangl2 to 
sites of cell-cell contacts [24]. We subsequently studied whether the 
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Figure 4. Effect of TG repeat variants on Celsr1-VANGL2 
interactions by Co-lp. Western blotting of the lysate with anti-GFP 
demonstrates the presence of wild-type GFP-Celsr1 protein and mutant 
proteins with different TG repeats. Following HA immunoprecipitation, 
blotting with anti-GFP confirmed physical association between Celsrl 
wild type and VANGL2, while this association was affected in CELSR1 TG 
repeat variants. 

doi:1 0.1 371 /journal.pone.0092207.g004 

C.5050-5051insTG and the C.5719-5720delTG affect the ability 
of CELSR1 to recruit VANGL2 to the cell-cell contact. MDCK II 
cells were transfected with PDS-RedCl-VANGL2 and GFP- 
CELSR1 wild type or its variants. CELSR1 and VANGL2 
subcellular localization were examined by confocal microscopy. As 
shown in Figure 3, wild type CELSR1 co-localized with VANGL2 
at the cell-cell contact in all of the cells examined (100%), while 
TG indelmutant forms of CELSR1 could not recruit VANGL2 to 
the contacting interface between co-transfected cells. 

CELSR1 C.5050-5051insTG and C.5719-5720delTG impair 
interaction with VANGL2 

We also tested whether the CELSR1 C. 5050-505 linsTG and 
the C.5719-5720delTG affect the interaction between CELSR1 
and VANGL2. HEK293T cells were co-transfected with HA- 
VANGL2 and either wild type or mutant GFP-CELSR1, 
incubated for 24 hrs and protein lysates were subjected to co- 
immunoprecipitation (Co-IP) and western blot analysis. The 
expression levels of mutant and wild type GFP-CELSR1 were 
found to be similar. Co-IP assay demonstrated that C.5050- 
5051insTG and C.57 19-5720delTG forms of CELSR1 impaired 
Celsrl interaction with HA-VANGL2 (Figure 4). 

Discussion 

Our study identified novel CELSR1 TG indels and SNV's in 
spina bifida patients. Two previous studies reported that rare 
mutations in CELSR1 are associated with human NTDs [19,20]. 
Both previous studies identified only SNVs. One study investigated 
the biological effects of some other SNV's in detail, so our study 



focused on evaluating the biological effect of the observed TG 
repeat mutations. 

Both the deletion and insertion were identified to be TG 
dinucleotide repeats. The insertion c. 5050-505 linsTG added a 
TG dinucleotide to TGTG, whereas c.5719-5720delTG removed 
a TG dinucleotide repeat from TGTGTG. It is known that among 
di-nucleotides, (TG)n are the most frequent in both humans and 
mice [25]. In humans, CELSR1 coding sequence region (CDS) has 
46 TGTG and 3 TGTGTG. In mice, Celsrl CDS region has 63 
TGTG and 6 TGTGTG repeats. In this California spina bifida 
cohort, the CELSR1 CDS TG dinucleotide repeat variation rate 
was 1% (2 in 192). 

Previous studies showed that missense mutations in membrane 
associated PCP genes including CELSR1 could affect PCP 
pathway signaling by disrupting membrane localization [19]. 
Here, micro insertions/deletions created truncated CELSR1 
forms lacking the transmembrane domain, so that membrane 
localization might be prevented. Indeed, absence of mutant 
CELSR1 from the membrane and its accumulation in the 
cytoplasm and nuclear compartments was observed. 

CELSR1 is known to physically associate with VANGL2, and 
CELSR1-VANGL2 interaction plays an important role in 
maintaining planar cell polarity during cell proliferation. During 
mitosis, CELSR1 interacts with VANGL2 and recruits VANGL2 
to endosomes. Following mitosis, PCP proteins are recycled to the 
cell surface, where asymmetry is re-established by a process reliant 
on neighboring PCP [16]. Mutations affecting CELSR1-VANGL2 
interactions can preclude VANGL2 recruitment to the endosomes 
during mitosis, thus disrupting PCP signaling. Both of the TG 
dinucleotide repeat variants identified in this study prevented 
Celsrl physical association with VANGL2. 

VANGL2 is another core PCP protein. It is a four-pass 
transmembrane protein and its proper localization is critical for 
VANGL2 function. Several proteins are required to establish 
VANGL2 localization. Previous studies demonstrated that NTD- 
inducing SNVs in Vangl2 itself, such as p.D255E and p.S464N, 
could lead to mislocation of the protein in mice [26,27]. SEC24B 
is a transport protein involved in vesicle trafficking, including 
shaping of the vesicle, cargo selection and concentration. 
Mutations of SEC24B affected VANGL2 membrane localization 
[28,29]. Here, we demonstrated that truncation of CELSR1 can 
disrupt VANGL2 cell-cell contact localization. Our results are 
consistent with the findings in the previous study by Devenport 
and Fuchs (2008) [24], which showed that deletion of the 
cytoplasmic tail of GFP-Celsrl partially impaired its localization 
to contacting interfaces and its ability to recruit Vangl2 at these 
cell-cell contacts. In our in vivo study, when GFP-Celsrl wildtype 
was transfected alone, Celsrl was distributed uniformly at 
membrane junctions (Figure 3A). When GFP-Celsrl was co- 
transfected with DsRedCl-VANGL2, it was distributed asymmet- 
rically to cell-cell contacts (Figure 3B). This observation is 
consistant with the Devenport and Fuchs (2008) study [24], which 
indicated that Vangl2 and Celsrl are dependent on one another 
for their proper asymmetric distribution. 

In summary, our study on a California spina bifida cohort 
indicates that mutations in CELSR1 contribute to the development 
of spina bifida. About 1% (2 in 192) of the spina bifida cases 
presented TG indels in CELSRL It is interesting that the two TG 
indels identified in this study caused severe biological malfunction 
of the mutant CELSR1 proteins, yet the birth defect associated 
with them is myelomeningocele. In contrast, the study by 
Robinson and co-workers [19] found SNV's that moderately 
altered the biological function of the CELSR1 protein, yet the 
associated defect was craniorachischisis, a more severe type of 
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NTD. One possible explanation is that NTDs are caused by a 
combination of multiple genetic and environment factors. One 
functional mutation is not sufficient to produce an NTD 
phenotype. There likely needs to be multiple functional mutations 
working together to cause an NTD phenotype, and to determine 
the NTD subtype. It was demonstrated in mice that homozygous 
PCP mutations caused craniochisichisis, as shown for Vangl2, 
Celsrl, and Pfk7 while double heterozygosity for a PCP mutation 
and non-PCP mutation can cause spina bifida, as shown for 
Vangl2/Dactl double heterozygotes. Perhaps the CELSR1 TG 
indels combined with other non-PCP mutations caused the 
myelomeningoceles that we observed, and that in the Robinson's 
study, the moderate CELSR1 mutations combined with other 
functional PCP mutations resulted in craniorachischisis. In our 
study, none of the CELSR1 TG indels were combined with 
mutations in other sequenced PCP genes including VANGL1, 
VANGL2, DISHEVELLED 1,DISHEVELLED2, DISHEV- 
ELLED3, FZD6, SCRIB and PTK7. Exome sequencing or whole 
genome sequencing approaches are warranted in order to examine 
the role of additional PCP genes such as PRICKLE, ANKRD6 
(ankyrin repeat domain 6; also known as DIVERSLN, the 
orthologue of diego), FUZ, and to detect non-PCP genes that 
may contribute to these NTD phenotypes. 

The four SNVs, which were predicted to be damaging by the 
four computational programs used in this study, were mapped to 
three different domains of CELSRl. Two of them, p.Alal023Gly 
and p.Ilell24Met were mapped to the eighth cadherin repeat 
domain of CELSRl, the same domain where two previously 
identified NTD mouse mutations (p.Aspl040Gly and 
p.Aspl 1 lOLys) were located. These two SNVs may affect the 
function of cadherin repeat domain which plays an important role 
in cell-cell contact. Variant p.Thrl362Met was predicted to 
change a threonine amino acid to a methionine amino acid. 
Thrl362 is a putative phosphorylation site of both protein kinase 
B (PKB) and G protein-coupled receptor kinase (GRK) based on 
analysis using KinasePhos 2.0 (http://kinasephos2.mbc.nctu.edu. 
tw/). Variant p.Arg2497Cys was mapped to the seven transmem- 
brane domain. The amino acid characteristics change from 
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